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^\The Antennae, NGC 4038/9, is the classic example of a pair 
of galaxies in gravitational interaction. Optical photographs 
^ ( e.g. Arp 1966; Schweizer 1978) show spectacular tails asso- 
I ciated with these galaxies, spanning a total angular diam- 
O eter of ~ 20', corresponding to 145 kpc at the distance of 
+3 25 Mpc (Tully 1988) assumed throughout this paper. The 
^ system is part of a modest concentration of galaxies called 
the NGC 4038 group, lying within the Crater Cloud, which 
itself lies about 7 Mpc below the Virgo Southern Extension 
(Tully & Fisher 1987). 

The Antennae has been the target of a sequence of dy- 
namical models, and as a result, its dynamical history is 
understood as well that of any other interacting system. 
Computer models such as Toomre and Toomre's (1972) sim- 
ple three-body model and Barnes's (1988) N-body model in 
which self-gravitating halo, disc and nuclear components are 
taken into account, have shown that the tails can be caused 
by the tidal interaction, and imply that they are approxi- 
mately 500 Myr old (Barnes 1988; Hut & Sussman 1987). 
The basic validity of these models has been confirmed by 
kinematical studies of both the neutral (van der Hulst 1979) 
and ionized (Amram et al. 1992) hydrogen gas. 

As one of the closest examples of an obviously inter- 
acting system, the Antennae has, over the years, been the 
target for high resolution multiwavelength studies of the ef- 
fects of tidal interaction. Optical Ha emission knots (Ru- 
bin, Ford & D'odorico 1970) and coincident powerful radio 
knots (Hummel & Van der Hulst 1986; Wolstencroft 1988) 



ABSTRACT 

We report the results of an X-ray spectral imaging observation of the Antennae with 
the ROSAT PSPC. 55% of the soft X-ray flux from the system is resolved into discrete 
sources, including components identified with the galactic nuclei and large Hll regions, 
whilst the remainder appears to be predominantly genuinely diffuse emission from gas 
at a temperature ~ 4 x 10 6 K. The morphology of the emission is unusual, combining a 
halo which envelopes the galactic discs, with what appears to be a distorted, but well- 
collimated bipolar outflow. We derive physical parameters for the hot gas in both diffuse 
components, which are of some interest, given that the Antennae probably represents 
an elliptical galaxy in the making. 
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are found throughout the central parts of both galaxies. Van 
der Hulst's (1979) 21cm line study showed that about 70% 
of the Hi in the system is associated with the luminous tails. 
Infrared images (Bushouse & Werner 1990) indicate that the 
nuclei of both galaxies are likely to be sites of very active star 
formation, and this is backed up by CO observations (Stan- 
ford et al. 1990), which show the presence of large concen- 
trations of molecular gas, not only at the two galactic nuclei, 
but also at the contact region between the two discs. 

The Einstein observation of the Antennae, (Fabbiano 
& Trinchieri 1983) was the major X-ray observation of an 
interacting galaxy system prior to ROSAT. The results were 
rather inconclusive, primarily due to the limited angular res- 
olution of the EinsteinYPC, however the emission was clearly 
extended, and contained a soft X-ray contribution not nor- 
mally seen in spiral galaxies, which the authors interpreted 
as probably arising from hot gas. Harder emission and a pos- 
sible hard point source in the vicinity of the contact region 
between the two galaxies were also seen. 

The ROSAT X-ray telescope, with the Position Sensi- 
tive Proportional Counter (PSPC) (Pfeffermann et al. 1986) 
in its focal plane, offers three important improvements over 
the Einstein IPC. It has much better spatial resolution, the 
90% enclosed energy radius at 1 keV being 27" (Hasinger et 
al. 1992). Since, as we will see below, the source is complex, 
and only a few arcminutes in diameter, this improvement 
is critical. Secondly, the spectral resolution of the PSPC 
(AE/E ~ 0.4 FWHM at 1 keV) is significantly better than 
previous imaging instruments, allowing characteristic source 
and diffuse emission temperatures to be derived. Thirdly, the 
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internal background of the PSPC is very low (~ 3 x 10 -5 ct 
s -1 arcmin -2 ; Snowden et al. 1994), which allows low sur- 
face brightness emission to be mapped. 



OBSERVATIONS, 
RESULTS 



DATA ANALYSIS AND 



NGC 4038/9 was observed with the ROSAT PSPC on four 
separate occasions for a total time of over 31000 s. Each 
of the four datasets was reduced separately using the Star- 
link ASTERIX X-ray analysis system and they were then 
combined, weighting each according to its exposure time to 
maximise the signal/noise. Although the spectral response 
of the PSPC varies with time, the variation across the four 
Antennae datasets was found to have a negligible effect on 
the fitted models, so we were able to use the response for 
the longest data block for the combined data. 

The data were first 'cleaned' by removing high back- 
ground periods - typically 2-3% of the data. They were then 
binned into a spectral image (which we will also refer to as 
the 'data cube') covering a 0.2° x 0.2° region and spanning 
the energy range 0.1-2.3 keV. A background model was con- 
structed using source-free regions of an annulus situated out- 
side the central PSPC support ring, plus knowledge of the 
energy dependent vignetting function, and was subtracted 
from the data cube. 

Individual ROSAT pointings are subject to an r.m.s. 
attitude uncertainty of approximately 6" (Hasinger et al. 
1992), hence before co-adding the four component spectral 
images, an attempt was made to coalign them by registering 
each on a nearby catalogued source (the F star HD104456). 
Since X-ray positional errors for this source were 6-9" (at 
90% confidence), the improvement in registration was only 
modest, and positions from the co-added data are subject 
to PS 4" r.m.s. systematic errors. 

From the co-added data cube, images were extracted 
in several different spectral bands. Fig. 1 shows contours 
of (0.1-1.9 keV) X-ray emission superimposed on an opti- 
cal (Schmidt telescope) image, whilst Fig. 2 shows similar 
images in three different energy bands. Each X-ray image 
has been lightly smoothed with a gaussian of standard devi- 
ation 10" to suppress noise. The contour levels in each plot 
increase by factors of two from 3.6 x 10 -4 ct s -1 arcmin -2 . 

The brightest X-ray emission covers the optical discs 
of the galaxies, with X-ray peaks centred on NCG 4038 (to 
the N) and NCG 4039 (to the S). These peaks are spectrally 
hard, whilst the remaining disc emission is more extensive in 
the two softer bands. A striking feature of Fig. 1 is the pres- 
ence of low surface brightness X-ray features extending to 
the north-west and south-west, and apparently culminating 
in point sources (labelled P and Q in Fig. 1). The emission 
from these X-ray extensions is mostly very soft. 

A map of the X-ray hardness is shown in Fig. 3. The 
grey-scale image indicates the average photon energy within 
each pixel containing a flux greater than 3.6 x 10 -4 ct s -1 
arcmin -2 (i.e. corresponding to the lowest contour levels of 
Fig. 1). The lightest regions correspond to a mean channel 
energy of ~ 0.2 keV, the darkest, to ~ 1.4keV. Whereas 
very hard emission is seen within the central parts of the 
southern disc, the central disc of NGC 4038 is rather softer 
(apart from a hard spot at its NE end). Much softer emission 




Right Ascension 

Figure 1. Contours of 0.1-1.9keV X-ray emission are shown su- 
perimposed on an optical image of NGC 4038/9. The contour 
levels increase by factors of two from 3.6 X 10 — 4 ct s — 1 arcmin -2 . 



predominates over the outer parts of both discs, and also 
in the two X-ray extensions, though small 'hot spots' are 
visible. 

A maximum entropy reconstruction (Gull 1989; Skilling 
1989) performed on a 0.1-1.9 keV image of the central regions 
of NGC 4038/9 and on soft (0.1-0.9 keV) and hard band (0.9- 
1.9 keV) images, is shown in Fig. 4. Seven discrete features 
(labelled A-G) can be distinguished. Best positions for these 
features were evaluated by locating peaks in the maximum 
entropy map (see Table 1). 

Firstly, an integrated spectrum for the system was ex- 
tracted from a circular region of radius 2.5', and fitted with 
standard spectral models (hot bremsstrahlung, power-law 
and Raymond and Smith (1977) hot plasma). The best fit 
was obtained with a Raymond and Smith hot plasma model 
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of temperature 0.69_ 03 keV and metallicity 0.03 
solar. The fitted column of 5. 5^ '\ x 10 20 cm -2 is signif- 
icantly greater than the column out of our own Galaxy 
(3.4 x 10 20 cm -2 ), indicating that either some fraction of the 
emission is intrinsically absorbed, or the emission is multi- 
component. The rather poor fit ( x 2 = 23.7 with 14 de grees of 
freedom) suggests that a separation of the emission into its 
constituent components (i.e. sources and diffuse emission) is 
necessary. 

For the purposes of further spectral analysis, the flux 
from point sources was removed and emission divided into 
two components - the bulk of the emission, which covers 
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Figure 2. Contours of X-ray emission in three different energy bands (a: 0.1-0.4keV, b: 0.4-0.9keV, c: 0.9-1.9keV) are shown super- 
imposed on optical images of NGC 4038/9. The contour levels in each of the plots increase by factors of two from 3.6 X 10 — 4 ct s — 1 
arcmin - 2 . 



the two discs, and the X-ray extensions to the north and 
south. A spectrum of the diffuse disc emission was obtained 
by extracting a circular region of radius 2.5' containing the 
disc emission from the background-subtracted data cube, 
and removing data at the positions of the seven features 
A-G. The remaining data were then collapsed into a spec- 
trum and exposure corrected. Spectra for the X-ray exten- 
sions were extracted in a similar way from approximately 
rectangular regions enclosing the lowest surface brightness 
contours shown in Fig. 1. In order to account for diffuse flux 
lost when removing the discrete sources, these spectra were 
renormalised using 'patched' images, in which the holes left 
at the source positions were filled by bilinear interpolation. 

The patched image referred to above is shown in Fig. 5. 
The seven features A-G have been removed, as have sources 
detected outside the inner regions of NGC 4038/9. Separa- 
tion of these features from the underlying emission proved 
difficult - a trade-off has to be made between removing all 
the source flux (implying patching to a large radius) and ob- 
taining a reliable estimate of the diffuse flux underneath it 
(which may require a smaller radius if the diffuse emission is 
structured). Features A-G were all removed to the 95% en- 
ergy enclosed radius for point sources. Sources outside the 
optical discs were treated similarly except in a few cases 
where a smaller radius was used to avoid biasing the resid- 
ual diffuse flux. The final map of the diffuse emission shows 
that it covers the whole of both optical discs, though not 
the tidal arms, and extends at least 2' (15 kpc) beyond the 
optical confines of the system to the NW and SW. Within 
the disc emission there appear to be excesses to the west of 
NGC 4038 and along the line where the two galaxies make 
contact. 

Extraction of spectra for the discrete sources A-G re- 
quired some care. The 95% enclosed energy radius of the 
PSPC on-axis point spread function (PSF) is PS 0.6' (at 
~ 0.9 keV), and it can be seen from Fig. 4 that circles of 



this size at the positions of the sources would overlap. Non- 
overlapping circles with radii given in Table 1 were centred 
on each source, and a spectrum extracted from each circle. 
A spectrum of the background plus diffuse disc emission was 
derived from the source free regions of the disc, and this was 
subtracted from each of the source spectra, with a normal- 
isation determined in each case by the flux underlying each 
source in the patched image. Finally, the resulting spectrum 
for each source was scaled up to allow for the flux lost in the 
wings of the PSF (see Table 1). This final scaling is really 
energy dependent, since the PSF width varies with photon 
energy, being larger at high and low energies. Since we ig- 
nored this effect we will have failed to scale up the lowest 
and highest energy emission sufficiently, relative to the cen- 
tre of the spectrum. In the worst cases, we calculate that 
we will have underestimated the flux by 50% at 0.2 keV, 
and by 20% at 2 keV. Simulations show that the low energy 
deficit would lead us to overestimate the column by up to 
20% and the high energy deficit would cause our fitted tem- 
perature to be 10-15% too low. These worst case distortions 
are comparable to the statistical errors quoted in Table 2. 

The spectra were fitted with standard spectral models: 
power law, thermal bremsstrahlung, and Raymond & Smith 
(1977) hot plasma models. Bremsstrahlung and hot plasma 
models gave results of similar quality, and the latter are 
generally quoted below, since they are physically better mo- 
tivated. Gaussian errors are not a very good approximation 
for these sources, so a maximum likelihood criterion, which 
allows explicitly for the Poissonian nature of the statistics, 
was used in preference to minimum chi-squared for the spec- 
tral fitting. 

The best fits to the source spectra are shown in Ta- 
ble 2. The comparative performance of different models can 
be judged on the basis of their likelihood ratios by the use of 
the Cash statistic (Cash 1979). In the case of the brightest 
sources A and D, for which the errors are fairly Gaussian, we 
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Table 1. Positions of sources A-G together with possible identifications. The radii used in the spectral extraction are given together 
with the net source counts within each of the circles, and the PSF correction factors used to scale up the spectra to allow for the flux 
lost in the PSF wings. 
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Figure 3. An X-ray hardness map of NGC 4038/9 superimposed 
on an optical outline. The average photon energy within each pixel 
containing a flux greater than 3.6 X 10 — 4 ct s — 1 arcmin - 2 (i.e. 
corresponding to the lowest contour level of Fig. 1) is shown. The 
lightest regions correspond to a mean channel energy of ~ 0.2 keV, 
the darkest, to ~ 1.4keV. 



were able to confirm that differences in the Cash statistic for 
different models corresponded closely to differences in chi- 
squared. All of the sources apart from E are better fitted by 
Raymond and Smith than by power law models. Where one 
fit is only marginally better than the other (Cash statistic 
values differing by < 2), the results for both are given, with 
the parameters of the poorer fit shown in brackets. 



The temperatures of sources A-D are very similar and 
are quite well constrained (to within ~0.25keV). Two lumi- 
nosities are shown for each source; one is the source luminos- 
ity which escapes from the galaxy, the other is the intrinsic 
source luminosity, i.e. before absorption in the host galaxy. 
For each object, we have also calculated a 'hardness ratio', 
equal to the difference divided by the sum of the high- (1.0- 
2.5 keV) and low- (0.1-1.0 keV) energy counts. This is useful 
for comparison with the work of other authors. The final 
column of the table lists possible counterparts for some of 
the sources, as discussed in the next section. 

The spectrum of the diffuse disc emission was also fitted 
with standard spectral models. The larger count rate in this 
case allowed the chi-squared statistic to be used. A Raymond 
and Smith model fitted significantly better than a simple 
power law spectrum (\ 2 = 21.6 with 14 degrees of freedom, 
compared to \ 2 = 35.7 with 15 degrees of freedom) and the 
results for this 'one component' model are shown in Table 3. 

The X-ray spectra of normal spiral galaxies are rather 
hard, and are believed to be dominated by emission from 
X-ray binaries (Fabbiano 1989). One expects that a similar 
contribution, due to unresolved hard sources in the bulges 
and discs of the two galaxies, should be present in the An- 
tennae. To investigate this, the diffuse spectrum was also 
fitted with a two component model, a Raymond and Smith 
hot plasma model (representing the truly diffuse emission) 
plus a hot (lOkeV) bremsstrahlung model (representing the 
contribution from unresolved X-ray binaries). The diffuse 
component column was constrained to be no less than the 
column out of our own Galaxy (3.4 x 10 20 cm -2 ). The best 
fit model parameters are shown in Table 3, and this model 
is overlaid on the diffuse spectrum in Fig. 6, whilst Fig. 7 
shows confidence contours for temperature against metallic- 
ity of the Raymond and Smith component. 

In Table 3, the total 0.1-2.0 keV luminosity of the dif- 
fuse disc emission is shown for both emission leaving the 
galaxy (escaping), and that emitted within the host galaxy 
(intrinsic). Only the hot bremsstrahlung component of the 
two-component fit is affected by the host column (its fitted 
column is significantly greater than the column out of our 
own Galaxy). 

The spectra extracted for each of the X-ray extensions 
were also fitted with standard spectral models, with the col- 
umn frozen at the Galactic value. Both northern and south- 
ern features were best fitted with hot plasma models, with 
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Table 2. Best fits to the spectra of discrete sources within NGC 4038/9, with one-sigma errors for one interesting parameter. Source 
spectra A, B, C, D, F and G are best fitted by Raymond and Smith hot plasma models. Source spectrum E is best fitted by a power-law 
model. Both escaping and intrinsic source luminosities are shown. The 'hardness ratio', is the difference divided by the sum of the high 
(l.0-2.5keV) and low (0.1-1.0keV) energy counts. Luminosities are based on an assumed source distance of 25 Mpc. 



Spectral parameters 


One-component 
model 


Two-component 
model 


Net counts (scaling factor) 


521±32 (1.78) 


Raymond Column (10 20 cm -2 ) 

and Smith Emission measure (10 60 cm -3 ) 

plasma Temperature (keV) 

(RS) Metallicity (Solar) 
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Hot Column (10 20 cm" 2 ) 

Brems. Emission measure (10 60 cm -3 ) 

(Br) Temperature (keV) 




87 5 + 215 
2200+ 4500 
10.0 (Frozen) 


X 2 (No. degrees of freedom) 


21.56 (14) 


17.80 (12) 


log-^X (escaping) 
log Lx (intrinsic) 


40.63 
40.63 


RS comp: 40.61 
Br Comp: 39.51 

Total : 40.64 
RS comp: 40.61 
Br Comp: 40.20 

Total : 40.75 



Table 3. Best fit parameters and luminosities resulting from fitting the diffuse disc emission spectrum. Also given is the number of 
counts extracted for the diffuse disc spectrum and the factor by which this is scaled to allow for the diffuse flux lost when removing the 
discrete sources (see text). Errors are one-sigma for one interesting parameter. 



very similar temperatures to the diffuse disc emission (see 
Table 4). 

Finally, we have investigated the spectral properties of 
the outlying sources P and Q. Power law models gave the 
best fits to both spectra, though thermal spectra are not 
ruled out. The best fit parameters are shown in Table 5. 



3 DISCUSSION 
3.1 The nuclei 

Millimetre wave CO observations show large quantities of 
molecular gas in the nuclei of both galaxies (Stanford et al. 
1990): 8xlO 8 M in NGC 4038, and 2xlO 8 M in NGC 4039. 
Optical and IR images indicate that both are sites of active 
star formation (Bushouse & Werner 1990), and star forma- 
tion rates of PS 2Mq yr -1 have been derived for both nuclei 



from corrected Ha fluxes (Stanford et al. 1990). They are 
also unusually bright in the radio - their extension and flat 
indices indicating a starburst origin for the emission, with 
inferred supernova rates of 0.14 yr -1 and 0.07 yr -1 respec- 
tively (Hummmel & Van der Hulst 1986). 

The X-ray feature B appears to be associated with the 
above-mentioned powerful radio, CO and Ha features ob- 
served at the centre of NGC 4038 (e.g. source 7 of Wolsten- 
croft 1988). However, it does show a significant offset - the 
centre of feature B lies 10" west of the (Stanford et al. 1990) 
CO feature and 13" west of the (Hummel & van der Hulst 
1986) radio feature. Similarly, feature D appears to be as- 
sociated with the powerful multi-wavelength sources at the 
centre of the southern galaxy (e.g. source 1 of Wolstencroft 
1988). Again offsets between the different features are ob- 
served, though they are somewhat smaller (feature D lies 9" 
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Table 4. Best fits to the northern and southern extension spectra with one-sigma errors for one interesting parameter. In both cases, 
Raymond and Smith hot plasma models gave the best fit. Luminosities are based on an assumed distance of 25Mpc. 
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Table 5. Best fit parameters for power-law fits to the spectra of sources P (to the north) and Q (to the south). Bremsstrahlung models 
gave marginally worse fits and the resulting temperatures are shown in brackets. Luminosities are based on an assumed source distance 
of 25 Mpc. Errors are one-sigma for one interesting parameter. 



north of the (Stanford et al. 1990) CO feature but only 5" 
north of the (Hummel & van der Hulst 1986) radio feature). 

The X-ray spectra of both 'nuclear' components are 
best fitted by absorbed, low temperature plasma models. 
The level of intrinsic absorption seen, however, is low com- 
pared to the optical extinction inferred for the nuclear Ha 
emission of ~ 2 — 3 magnitudes, equivalent to a column 
Nh ~ 5 x 10 21 cm -2 . This, together with the ~ 1 kpc offsets 
in X-ray emission compared to CO and radio emission, sug- 
gests that we are not seeing the nuclei themselves in X-rays 
but instead are observing hot gas near the nuclear regions. 
Such extranuclear hot gas might be expected given the large 
star-formation rates and extended flat index radio features 
(indicative of a great deal of supernova activity) observed in 
the nuclei. The fact that the southern nucleus has the higher 
soft X-ray luminosity, whilst the northern one contains more 
molecular gas and has a higher star formation rate and ra- 
dio flux, also supports the idea that we are not seeing the 
central starbursts directly. 

The soft spectra of these features (hardness ratio ~ -0.2 
to -0.3) are comparable to those of superbubbles seen in the 
LMC (Williams & Chu 1995), though their luminosities are 
two orders of magnitude higher. 



3.2 Emission from Hn regions 

The star-forming knots seen in Ha (Rubin et al. 1970) are 
predominantly found in NGC 4038 and correspond to mas- 
sive (~ 10 7 Mq) Hll regions. The total star-formation rate 
over all the knots can be estimated by combining the calcu- 
lated nuclear star-formation rates of Stanford et al. (1990) 
with the total Ha luminosity one expects to observe based on 
the thermal radio emission (Hummel & van der Hulst 1986). 
This global star-formation rate comes to ~ 20MQyr -1 , an 
order of magnitude more than the total in our own Galaxy. 
Corresponding knots are seen in the radio (see Wolstencroft 
1988) and contain about 35% of the total 1.5 Ghz emission 



(Hummel & van der Hulst 1986). Some of these radio knots 
are actually brighter than those at the galactic cores. 

Both X-ray features A and C seem coincident with op- 
tical and radio counterparts which are most readily inter- 
preted as large Hll regions. The centre of source A, the 
brightest of the X-ray features, lies within the largest of the 
Ha knots (source E of Rubin et al. 1970), and ~8" north- 
east of the centre of radio knot 5 (Hummel & van der Hulst 
1986). Source C lies very close (within 2-3") to radio knot 
10 (Hummel & Van der Hulst 1986) and Ha knot T (Ru- 
bin et al. 1970). The radio sources associated with A and 
C have the flattest radio indices of all the observed radio 
knots, implying the highest ratio of thermal to non-thermal 
emission (Hummel & van der Hulst 1986). This suggests 
strongly that the X-ray emission from A and C arises from 
collections of supernova remnants. This idea is further sup- 
ported by the fact that a comparison of the radial profiles 
of features A and C with the Point Spread Function (PSF) 
of the ROSAT PSPC indicates that both features appear 
marginally extended on a scale of between 1.5 and 4.5 kpc. 
However, given the presence of underlying diffuse emission, 
higher resolution X-ray observations (which should be avail- 
able soon from the ROSAT High Resolution Imager) are 
required to confirm this. 

Again, as with features B and D, the hardness ratios 
of these features can be compared with ratios of known fea- 
tures in the LMC and in M101 (Williams & Chu 1995). As 
with features B and D, the ratios are low (though not as 
low), and are comparable with features such as the old, ex- 
tended supernova remnant, N44, in the LMC. The luminosi- 
ties of features A and C are higher, by a factor of 3-30 than 
the brightest features seen in M101, and, if we have truly 
resolved features A and C, they are also larger (the Hll re- 
gions within M101 have diameters up to 1 kpc (Kennicutt 
1984)). 

Feature A may well be associated with the near-IR 
bar-like feature extending from the northern nucleus to the 
south-east to radio knot 5 (close to X-ray source A) (Wright 
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Right Ascension 

Figure 4. Maximum entropy reconstructions of NGC 4038/9 in 
three different energy bands (a: 0.1-1.9keV, b: 0.1-0.9keV, c: 0.9- 
1.9keV) superimposed on optical images. The contour levels in- 
crease by factors of two. 
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Figure 5. Residual (0.1-1.9keV) diffuse emission from NGC 
4038/9 after removal of sources, and linear interpolation over the 
resulting 'holes'. Contour levels increase by factors of two from 
3.6 X 10 _4 ct s _1 arcmin -2 . 



& Mclean 1987). It is suggested that the non-axisymmetric 
gravitational potential associated with this bar might be re- 
sponsible for driving gaseous material towards the nucleus. 
It is now well established that substantial starburst emission 
can occur at the ends of bars as well as at the nucleus (e.g. 
Knapen & Beckman 1994). 

3.3 Other discrete X-ray features 

Sources E, and F, lying south-west of the southern nucleus, 
and source G, lying approximately due north of the northern 
nucleus, are all apparently, and intrinsically, less luminous 
than sources A-D discussed above. Source E is the only one 
of these with a possible counterpart at other wavelengths - 
it may be associated with a weak optical knot (source AA) 
noted by Rubin et al. (1970). 

Sources E and G lie within the optical confines of the 
galaxies, and the fits to the spectra of both indicate absorb- 
ing columns ~ 10 21 cm -2 ), well above the column out of our 
Galaxy. The most likely candidates for foreground Galactic 
sources at these flux levels, late type stars, can therefore be 
ruled out, since their X-ray spectra do not show intrinsic 
absorption. It is also very unlikely than E and G could be 
background objects (e.g. quasars) shining through the entire 
Hi column of the Antennae. Using the quasar log N — log S 
function of Hasinger et al. (1993), the probability of even one 
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Figure 6. Diffuse emission spectrum from NGC 4038/9 with best 
fit two-component model (solid line). The contributions of the 
two components are shown as dashed (Raymond and Smith) and 
dotted (hard bremsstrahlung component) lines. 

Chi— squared 

I 1 1 1 1 1 




0.2 0.4 0.6 

Temperature (keV) 



Figure 7. 99%, 95%, 90% and 68% confidence contours in the 
temperature-metallicity plane for the Raymond and Smith com- 
ponent of the two-component fit to the diffuse disc spectrum (see 
Table 3). All other parameters (except the hot bremsstrahlung 
temperature, which is frozen at lOkeV) are left free to optimise. 



background quasar appearing within the optical confines of 
the system at flux levels as high as those observed, through 
an absorbing column consistent with the fitted columns of 
sources E and G is only 0.02. 

It seems likely then, that E and G are sources within 
the Antennae. Source E is actually the hardest of all the 
point sources, and stands out clearly in Fig. 3. It is not as 
hard as known X-ray binaries within the LMC, e.g. LMC 
XI and R136 (Williams & Chu 1995), but is harder than 
the superbubble features discussed earlier, and may there- 
fore be a collection of hot supernova remnants with possibly 
some contribution from X-ray binaries. Approximately 50- 
100 sources similar to LMC X-l would be needed to explain 
the high X-ray luminosity of source E. Source G is likely, 
given its low fitted temperature, to be a superbubble or a 
collection of supernova remnants, though its X-ray luminos- 
ity is some two orders of magnitude greater than superbub- 
bles in the LMC (Williams and Chu 1995). 

Source F lies outside the optical confines of the system 
and has no optical counterpart (the nearest optical source 
with a magnitude < 20 m is almost 0.5' away). The X-ray 
source is weak, and its spectrum is poorly constrained, but 
clearly is very soft. Bright 'supersoft' sources, which may 
have luminosities in excess of 10 38 erg s -1 have been de- 
tected by ROSAT in a number of nearby galaxies (e.g. Ka- 
habka, Pietsch & Hasinger 1994), but in the present case 
the fact that the source lies just outside the visible disc of 
NGC 4039 makes such an identification unlikely. The soft 
spectrum is consistent with a foreground white dwarf, but 
this would be expected to be brighter than ~ 17 m in the 
optical. Hence, assuming that the feature is a point source, 
it seems most likely (given that it lies beyond the absorbing 
effects of the galactic discs) to be an unusually soft, optically 
faint background quasar. 

3.4 Residual disc emission 

The residual emission covering the optical discs, after re- 
moval of features A-G, constitutes ~40% of the total flux 
and appears to envelope the whole system apart from the 
tidal tails. Some structure in this emission can be seen in 
Fig. 5. The excess emission on the western side of NGC 4038 
follows closely the arc-like structure of knots seen in both 
the radio (Hummel & van der Hulst 1986) and Ha (Rubin 
et al. 1970), and peaks close to Ha knot G (Rubin et al. 
1970) and radio knot 13 (Hummel & van der Hulst 1986). 
In addition, there appears to be enhanced emission stretch- 
ing east-west along the interface between the two galactic 
discs and peaking at the eastern end. 

The best one-component fit to the residual disc spec- 
trum is a low temperature (0.49 keV) Raymond and Smith 
plasma, absorbed by a column approximately equal to the 
Galactic hydrogen column (3.4xl0 20 cm -2 ). Adding an ad- 
ditional 10 keV bremsstrahlung component to represent un- 
resolved hard sources improves the fit somewhat (Table 3). 
However, as can be seen from Fig. 6, this hard component 
only contributes significantly towards the top of the PSPC 
energy range, and an F-test shows that the improvement in 
fit quality is significant at less than 90% confidence. Hence 
although such a hard contribution is to be expected, it is 
not clear that we have detected it. The error region for the 
parameters of the hard component allowed by our data are 
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correspondingly large (Table 3) - a higher luminosity is al- 
lowed provided that the intrinsic absorption associated with 
the component is increased, such that it continues to have 
little impact in the ROSAT energy range. For the best fit 
model, the hard component contributes ~8% of the 'dif- 
fuse' disc flux, and has an intrinsic 0.f-2.0keV luminosity of 
f.6 x fO 40 erg s -1 , and a bolometric luminosity ~ 3.5 times 
larger. This compares, for example, to a total luminosity in 
discrete X-ray sources of ~ 2.6 x f 39 erg s -1 in our own 
Galaxy (Watson 1990). 

The luminosity of the soft diffuse disc component 
(which is affected very little by uncertainties over the hard 
component) is 4.1 x 10 40 erg s _1 in the 0.1-2.0 keV band. 

3.5 The collision interface 

Direct evidence for interaction of the two galaxies in seen 
in the Antennae in the form of a large molecular gas com- 
plex (seen in CO emission) and a peak in the diffuse radio 
flux (Stanford et al. 1990; Hummel & van der Hulst 1986) 
at the interface between the discs. The high magnetic field 
(ss 40 fiG) implied by the radio observations probably arises 
from compression of the ISM. The enhancement in the resid- 
ual X-ray emission between the two discs seen in Fig. 5 has a 
maximum at the eastern end, close to the CO and diffuse ra- 
dio peaks, and is therefore also probably associated with the 
collision. This is not the same as the much more luminous 
hard source tentatively identified with the contact region by 
Fabbiano & Trinchieri (1983) on the basis of an observa- 
tion with the Einstein IPC. The higher spatial resolution 
of R OS A T shows that this hard emission actually originates 
predominantly from the powerful, absorbed source D, at the 
nucleus of NGC 4039. 

The hardness map, Fig. 3, shows that the emission from 
the interface region is harder than the bulk of the diffuse 
disc emission. An X-ray spectrum was extracted from an 
area of 0.6 arcmin 2 centred on the 'contact peak', and fit- 
ted with a hot plasma model, freezing the column at the 
value obtained from fitting the whole of the diffuse disc 
emission (3.3 x 10 20 cm -2 ). This resulted in a temperature of 
0.943.Q Q8 keV, significantly hotter than the remaining diffuse 
emission. If this is interpreted as collisionally heated gas, the 
temperature implies a shock speed of ~ 700 km s" 1 . Fig. 3d 
of Barnes (1988) indicates that the two galaxies should have 
a radial velocity difference of 300-400 km s -1 . Assuming a 
similar tangential component (not unreasonable, given the 
modelled orientation and history of the system), this trans- 
lates into a collision velocity of 450-550 km s -1 , in good 
agreement with inferred collision velocities for 'typical' in- 
teracting systems (~500 km s _1 ; Harwit & Fuller 1988). Al- 
ternatively, this 'hot spot' might be due to enhanced star 
formation taking place in the dense molecular gas near the 
contact region. 

3.6 The X-ray extensions 

The extensions of the diffuse emission to the NW and SW 
of the system which can be clearly seen in Figs.l and 5, look 
rather like X-ray analogues to the tidal arms. However they 
apparently have no counterparts at other wavelengths. This 
implies that they cannot have been separated from the main 



body of the system by gravitational forces, and cannot arise 
from unresolved discrete sources. 

The idea that we are seeing streamers of hot gas is sup- 
ported by the spectral properties of the X-ray emission (see 
Table 4). This hot gas has presumably been ejected from 
the body of the system by gas pressure forces. Two possible 
sources of energy for this are potentially available - injec- 
tion from stars and supernovae, and the kinetic energy of 
the galaxy-galaxy encounter. 

Starburst activity such as is seen in the Antennae is 
known in many cases to trigger galactic winds: massive 
galactic-scale outflows of hot gas driven by supernovae and, 
to a lesser extent, stellar winds (Heckman, Lehnert & Armus 
1993). The extended soft X-ray emission seen above the discs 
of many starburst galaxies both with Einstein (e.g. Wat- 
son, Stanger & Griffiths 1984; Fabbiano 1988) and ROSAT 
(Heckman 1993; Read 1994) are thought to be due to am- 
bient interstellar gas clouds in the galactic disc and halo, 
which have been shock-heated by a hot, tenuous wind. 

The temperatures and luminosities inferred for the X- 
ray extensions in the Antennae are similar to what is seen 
in the winds of the nearby starburst galaxies NGC 253 and 
NGC 3628 (Read 1994). However, the soft X-ray emission 
is visible to a larger radius (approximately double) in the 
Antennae. This may reflect the fact that the hot wind in 
the Antennae has encountered, and shock heated, cool gas 
which had already been ejected from the galactic discs as a 
result of the violent encounter. The presence of such ejecta is 
supported by the observation of a considerable amount of Hi 
beyond the optical confines of the galaxies (van der Hulst 
1979), and tidally stripped material apart from the main 
tidal tails is also seen in dynamical models of the system 
(Hut & Sussman 1987; Barnes 1988). A hot galactic wind, 
travelling at a velocity of ~ 2000 — 3000 km s -1 , can easily 
reach and shock high altitude clouds within the timescale 
(~ 10 8 yr) since pericentre indicated by modelling (Barnes 
1988). 

The morphology of the X-ray extensions in the Anten- 
nae differs from that of the classic bipolar galactic winds 
which emanate perpendicularly to the galactic disc in sys- 
tems such as M82 and NGC 253. Although the extension 
consists of two branches, these are not oppositely directed, 
but at an angle of ~ 120°. However the Antennae is a very 
rapidly evolving system. Over the past 2 — 3 x 10 8 years it 
is believed that NGC 4039 has swung round from the north 
and passed in front of NGC 4038 to reach its current posi- 
tion to its south. Substantial slewing of the inner galactic 
discs of both galaxies has also occurred. It is not too sur- 
prising then, if the wind departs from the standard bipolar 
form. We cannot even be certain, given the fact that there 
are two active star-forming nuclei, that we are not seeing 
two separate winds. 

The other possibility for powering a hot gas outflow 
is the direct thermalisation of some of the kinetic energy 
of the two colliding galaxies. Harwit et al. (1987) have in- 
voked direct collisions between molecular clouds to explain 
the very high infrared luminosities of some merging galaxies. 
On a larger scale, hydrodynamical models of merging galaxy 
clusters have shown (Evrard 1990) that gas is shock heated 
along the interface between the intracluster gas associated 
with each cluster, and that this hot, overpressured gas can 
then squirt out perpendicular to the collision axis. In order 
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for such a phenomenon to explain the rather linear feature 
seen in the Antennae, we would have to be viewing the disc 
of squeezed gas approximately edge-on. This conflicts with 
the picture which emerges from the dynamical modelling of 
the system (Barnes 1988), which implies that the angle be- 
tween the viewing direction and the collision axis is quite 
small (~ 30°). 

3.7 A pair of satellite sources 

One striking aspect of the X-ray extensions seen in Fig. 1 
is the fact that both appear to culminate in bright sources. 
These sources, which we designate P (to the NW) and Q (to 
the SW) are symmetrically disposed at a distance of ~20 kpc 
from the optical boundary of the galaxies. Both have ra- 
dial profiles consistent with the point spread function of the 
PSPC and have no catalogued counterparts. Distortion of 
the X-ray contours due to these sources can be seen in the 
Einstein image (see the figure in Fabbiano, Kim & Trinchieri 
1992), though they have never been commented upon. 

We have obtained an optical image from the photo- 
graphic co-addition of 5 sky-limited IllaJ Schmidt plates 
of the Antennae field to search for optical counterparts to 
the two sources. A stellar image with Bwl9.8 lies within the 
90% confidence error circle of P, and can be seen in Fig. 1, 
however nothing is visible within the error circle of Q down 
to B«24. 

Source P has an unusual spectrum, lacking in both high 
and low energy flux, which is fitted by a steep spectrum with 
an absorbing column which exceeds that out of our Galaxy 
(see Table 5). Source Q has a flatter spectrum and smaller 
column. We derived fluxes for the two sources at each of the 
four epochs at which our data were collected. The result- 
ing light curves are shown in Fig. 8. The errors plotted are 
purely statistical, but since the sources are not significantly 
contaminated by diffuse emission we expect systematic ef- 
fects to be negligible. Fitting the light curves to constant 
flux values results in a reduced % 2 of > 4 for both P and Q. 
Both sources appear to have varied in brightness by a factor 
of a few on a timescale of ~ 1 yr. This confirms that they 
must be point-like. 

A rather similar situation has recently been noted in 
NGC 4258 by Pietsch et al. (1994), who detect two bright 
pointlike X-ray sources, located symmetrically at a distance 
of 9 kpc from the weak Seyfert nucleus, and coaligned with 
the tips of a pair of 'anomalous arms' within the galactic 
disc. A mechanism for ejecting massive (~ 1O 8 M ) objects 
from galactic nuclei has been suggested by Saslaw, Valtonen 
& Aarseth (1974). If P and Q represented ejected objects, 
and if the ejection process could also give rise to a trail of 
hot gas, this would explain the curious juxtaposition of the 
sources and the X-ray extensions. 

Unfortunately, it is difficult to make this interesting idea 
work, given the substantial width (~ 7 kpc) of the exten- 
sions. A very high velocity projectile might generate a broad 
feature after adiabatic expansion from a narrow hot wake, 
but no temperature gradient is seen along the extensions (see 
Fig. 3). Alternatively, ejected objects might create a channel 
out of the system through which hot gas generated by star- 
burst activity could escape as a wind. However, the fact that 
the point sources appear at the ends of the X-ray extensions 
has no natural explanation within such a model. 




Figure 8. Variation in count rate from sources P and Q over the 
four separate observations. Errors are 1<t. 



Despite their suggestive symmetry and positioning rela- 
tive to the X-ray 'arms', P and Q might be entirely unrelated 
background or foreground sources. In this case what could 
they be? The most common sources in ROSAT fields at 
these flux levels are quasars. The mean spectrum of quasars 
in the ROSAT band can be represented by a power law 
with a photon index of 2.2 ± 0.2 (Branduardi-Raymont et 
al. 1994; Roche et al. 1995), though some individual spec- 
tra are known to be shallower (a ~ 1.6 — 1.7; Biihler et 
al. 1995). This is consistent with the spectral properties of 
source Q (its best fit column is low, but is consistent with 
the galactic value within 90% confidence), though it would 
be an unusually optically faint example (L.R.Jones, private 
communication). The X-ray spectrum of P is steeper than 
that of quasars. Its fitted temperature (0.37 keV) and possi- 
ble optical counterpart are consistent with a late-type star, 
but the inferred column is too large. Optical spectroscopy 
of both sources is required to make further progress. 



3.8 A tidal dwarf galaxy 

Mirabel, Dottori & Lutz (1992), have discovered a dwarf 
galaxy in the process of forming near the tip of the southern 
tidal tail of the Antennae. The nascent galaxy consists of 
a chain of nebulae ionized by very young (< 2 — 6 Myr) 
massive stars, embedded in an envelope of Hi gas and low 
surface brightness optical emission. 

This region lies outside the field shown in Fig. 1, but 
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well within the PSPC field. We detected no significant X- 
ray emission from this region, and derived a 2a upper limit 
to the count rate within a circle of radius 0.6', (the 95% en- 
closed energy radius of the PSPC point spread function for 
a 0.5 keV source) at the position of the massive (~ 10 9 Mq) 
Hi complex at the tip of the tidal tail (van der Hulst 1979). 
Assuming a 5 keV bremsstrahlung model and Galactic ab- 
sorption, this upper limit corresponds to a 0.1-2.0 keV lu- 
minosity Lx < 4.3 x 10 38 erg s -1 (the equivalent 1 keV 
bremsstrahlung upper limit is 4.0 x 10 38 erg s -1 ). Including 
an additional intrinsic column, due to the local Hi would 
raise these limits only slightly, since the column through the 
Hi complex is only 1-2 x 10 20 cm -2 . 

Dwarf irregular galaxies have a wide range of X-ray lu- 
minosities, but the fact that we detect no X-ray emission 
in this case, despite the active star formation, is not sur- 
prising if all the stars have been formed within the last 
6 Myr (Mirabel et al. 1992). Even the most massive of these 
stars will not have had time to evolve off the main se- 
quence into the X-ray bright supernova remnant or high- 
mass X-ray binary phase. From Mirabel et al. (1992), one 
can estimate a visual luminosity Lv for the dwarf galaxy of 
PS 3.8 x 10 41 erg s -1 . This gives to rise to an upper limit to 
Lx I Lv of 1.1 x 10 -3 . O stars have a maximum Lx /-^bol ra ~ 
tio of only 10" 5 44 (Chlebowski, Harnden & Sciortino 1989), 
and so would not be detectable in our X-ray data. 

3.9 Physical properties and history of the hot gas 

Mean physical properties for the hot gas can be inferred 
from our spectral fits if we make some assumptions about 
the geometry of the emission. Considering first the diffuse 
emission covering the discs of the two galaxies, this occu- 
pies a roughly circular region on the sky, and the simplest 
assumption to make is that it in fact constitutes a spherical 
bubble of the radius observed (~15kpc). This is the largest 
volume over which the gas is likely to be distributed. A ra- 
dial profile of the X-ray surface brightness, after removal of 
the discrete sources (i.e. profiling the image in Fig. 5), shown 
in Fig. 9 casts some doubt on the validity of the spherical 
bubble model. The surface brightness falls rapidly outside 
the region covered by the optical discs, suggesting that the 
diffuse emission may be bound to the individual discs in the 
form of hot coronae. Assuming an isothermal gas at the ob- 
served temperature (4 x 10 6 K)) and adopting an analytical 
Plummer-Kuzmin form (Binney & Tremaine 1987) for the 
potential of each disc, one derives a scale height of approxi- 
mately 4.8 kpc for the gas above each disc. This agrees quite 
well with the scale on which the surface brightness drops off 
outside the disc region in Fig. 9. Taking this as the depth 
of the gas layer along the line of sight, and assuming a sim- 
ple cylindrical geometry gives the minimum volume for the 
emitting gas. 

Using the volumes derived for the bubble and coronal 
models, the fitted emission measure r]n 2 e V (where r\ is the 
'filling factor' - the fraction of the total volume V which is 
occupied by the emitting gas) can be used to infer the mean 
electron density, n e , and hence the total mass Mgas, thermal 
energy and cooling time t coo i of the gas, and also the 
mass cooling rate M coo i and adiabatic expansion timescale 
^exp- 

Similar estimates were made (assuming uniform density 



and cylindrical geometry in the plane of the sky) for the 
X-ray extensions to the north and south. All the resulting 
parameters are listed in Table 6. 

In the case of the X-ray extensions, the morphology 
strongly suggests that there is some sort of collimated flow. 
If r] ~ 1 for the X-ray emitting gas, then they are presum- 
ably streamers of gas at T ~ 4 x 10 6 K flowing away from the 
system. Unless some additional force is confining them, these 
streamers will expand at their sound speed (~ 300 km s -1 ) 
and dissipate on a timescale ~ 3 x 10 r yr unless they are 
continually replenished. Alternatively we might be seeing 
soft X-ray emission from clouds with r] <C 1 which are be- 
ing shock heated by a fast collimated wind. Theoretical 
modelling of starburst winds tends to support this option 
(Suchkov et al. 1994), since it suggests that the soft X-ray 
luminosity of such clouds should considerably exceed that 
of the hot wind fluid itself. 

Turning to the gas which covers the discs, if this is even 
approximately in hydrostatic equilibrium, then its filling fac- 
tor is unlikely to be much less than unity, since it would then 
require hotter low density gas to fill the rest of the volume, 
and this could not be bound to the galaxies (which have 
typical escape temperatures of ~ 10 r K). The alternative is 
again that we are seeing gas clouds which are shock heated 
by an escaping hot wind, invisible in our waveband, in which 
event r] <C 1. Since the diffuse disc emission appears to en- 
velope the whole system, this would require either that the 
wind is directed quite precisely towards us, or that it is es- 
sentially spherical. The former case seems unlikely given the 
far from face-on aspect of the galactic discs, whilst for a 
spherical wind it is hard to understand the abrupt cutoff in 
emission at a radius of ~ 15 kpc, except in the direction of 
the two extensions. This cannot be due to a lack of material 
in the vicinity available for shocking by the wind, since for 
example there is plenty in the nearby tidal tails. Hence one 
would have to suppose that the hot spherical wind does not 
extend beyond the X-ray bright region. Since its velocity 
is expected to be ~ 2000 — 3000 km s -1 , this implies that 
it must have been initiated within the past 10 r yr, which 
seems unlikely since the interaction timescale is an order of 
magnitude longer. 

Assuming, as seems probable, that the diffuse disc gas 
is bound to the system and has r] ~ 1, its total mass, 
~ 1 — 1.5 x 10 9 Mq, is substantial. For comparison, the main 
disc region contains ~ 2 x 10 9 Mq of molecular gas (Stanford 
et al. 1990) and ~ 1 x 10 9 M Q of Hi (van der Hulst 1979). 
Assuming that this gaseous halo has been generated over the 
past 1 — 2 x 10 8 yr, the gas loss rate is ~ 8Mq yr -1 whilst 
cooling is currently returning the gas to the central regions 
at only a fraction of this rate. The energy requirement in- 
volved in generating this halo of hot gas is essentially equal 
to its thermal energy content (the gravitational potential en- 
ergy adds a further ~10%), which corresponds to the energy 
released in a few million supernovae. These figures may be 
compared with estimates for the total star formation rate 
in the system of 20Mq yr -1 , and for the supernova rate of 
~ 1.4 yr -1 , about 15% of which takes place in the nuclei 
(the supernova rate can be calculated using the rates given 
in van den Bergh & Tammann (1991) and scaling these with 
the far-infrared luminosity as in Devereux & Eales (1989)). 
Hence, mass and energy requirements of producing such a 
halo do not present a problem. 
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Total diffuse disc emission Northern Southern 



(bubbfe) (coronaf) efongation efongation 



log-^X (erg s 1 ) 


40.61 


39.30 


39.33 


Tx (keV) 
n e (cm -3 ) 


0.36 


0.44 


0.43 


0.0049 X 1/y/rj 


0.0077 X 1/y/rj 

9.7 X io 8 y^ 


0.0066 X 1/y/rj 


0.0042 X 1/y/rj 


M gas (M ) 


1.5 X w 9 ^/v 


7.1 X 10 7 ^ 


9.0 X W 7 ,fn 


^th ( er g) 


3.1 X 10 57 y5J 


2.0 X 10 57 ^J 


1.7 X 10 56 y5J 


2.2 X 10 56 ^J 


*cool (Myr) 


2410^ 


1530^ 


2700^ 


3260^ 


M cooY (M yr" 1 ) 


0.62 


0.63 


0.03 


0.03 


texp (Myr) 


59 


16 


29 


30 



Table 6. Values of physical parameters for the diffuse disc gas (using both hemispherical bubble and coronal models), and for the 
northern and southern extensions (see text), rj is the filling factor of the gas. 



The cooling time of the halo is over 10 yr, which sug- 
gests that old merger remnants might retain an X-ray bright 
halo for some time after the optical merger (believed in 
a merger such as the Antennae's to result in an ellipti- 
cal galaxy (Toomre 1977)) is complete. Such halos are, of 
course, commonly observed in elliptical galaxies (Fabbiano 
1989), but are attributed to material ejected during the early 
star forming history of early type galaxies. The halo we 
see in the Antennae differs significantly from those seen in 
elliptical galaxies. Fig. 9 shows that the diffuse X-ray flux 
falls sharply outside the disc region, with surface brightness 
S(r) oc r -4 ' 5 , and is contained within a radius r ~25kpc, 
apart from the X-ray extensions. In contrast, elliptical halos 
have X-ray surface brightness profiles which are much flatter 
(typically S(r) oc r -2 ), extend over radii ~ 50 — 100 kpc, and 
contain ~ 10 9 - 10 10 M Q of gas (Fabbiano 1989; Forman et 
al. 1993; Trinchieri et al. 1994; Rangarajan et al. 1995). The 
temperature of the hot halos in ellipticals is T fa 0.7—1.2 keV 
(Matsushita et al. 1994; Rangarajan et al. 1995), consider- 
ably higher than the T fa 0.4 keV seen in the Antennae. 

These differences do not preclude the possibility that 
the corona seen in the Antennae might evolve into a more 
massive, hotter and extended halo if more energy and mass 
is injected into it as the merger proceeds, over the next few 
xlO 8 yr. As is well known from the study of galaxy clus- 
ters (e.g. Gorenstein et al. 1978), an isothermal halo of hot 
gas in hydrostatic equilibrium in a potential well has a den- 
sity profile, and hence an X-ray surface brightness profile, 
determined by the ratio of its temperature to the virial tem- 
perature in the potential. The low temperature of the corona 
at present therefore leads directly to a steep profile, and the 
injection of extra specific energy, raising the temperature to 
values similar to those seen in ellipticals, would result in a 
much flatter profile (doubling the temperature would cause 
the surface brightness index to drop by almost a factor of 
two). 

It is also interesting to compare our results with the 
properties of the X-ray faint, early-type galaxy NGC 4365 
(Fabbiano, Kim & Trinchieri 1994; Pellegrini & Fabbiano 
1994), which is believed, on account of its counter-rotating 
core, to be the product of a recent merger. A very soft 
(T fa 0.2 keV) component is observed in the X-ray spec- 
trum of the system, accounting for about half of the (0.1- 
2.0 keV) flux. There is tentative evidence that the X-ray sur- 
face brightness extends beyond the optical profile, to a ra- 
dius of ~ 20 kpc, and Pelligrini & Fabbiano (1994) suggest 



that the observed soft spectral component arises from a com- 
bination of stellar emission from the galaxy, and a gaseous 
halo with T fa 0.6 keV, and L x ~ 2 x 10 40 erg s _1 . If this 
interpretation is correct, then the temperature and extent 
of this halo are comparable to what we see in the Antennae. 
Its luminosity is somewhat lower, but NGC 4365 is also an 
optically fainter system. 

The mass and thermal energy content of the X-ray ex- 
tensions is considerably smaller than that of the halo gas 
even if r] ~ 1 for this gas. However this is probably gas 
which is being lost entirely from the system. If we are seeing 
the wind itself (i.e. r] = 1) then the mass loss rate can be 
calculated from the inferred gas density, assuming a wind ve- 
locity. This gives a mass loss rate, including both extensions, 
of M ~ iooo Mq yr -1 , where «iooo is the velocity in 

units of 1000 km s - , which could significantly deplete the 
interstellar gas content of the galaxies if it continued for 
~ 10 8 yr. Similarly, the rate at which kinetic energy is lost 
in the wind is i w ind ~ 5 x 10 42 T/»??»iooo er g s_1 (both these 
results are maxima in the case where the extensions lie per- 
pendicular to our line of sight). Hence a significant fraction 
of the total supernova luminosity of ~ 4 x 10 43 erg s -1 may 
be carried off in the wind. If the hot wind itself is invisible to 
ROSAT then the above results still apply for the mass and 
energy lost in the clouds (where now we will have r] <C 1), 
provided that they do escape from the system. However, 
there is additional mass and energy loss associated with the 
hotter gas. 

As is shown in Table 3, the metallicity of the hot disc 
gas inferred from hot plasma fits is very low, and the results 
of fits to the X-ray extensions given in section 3.6 show the 
same phenomenon. This result seems very surprising given 
the high supernova rate inferred in the Antennae and the 
consequent injection of metal-rich gas expected. It is neces- 
sary to be cautious in interpreting metallicity results from 
ROSAT spectra, since the low spectral resolution makes it 
impossible to distinguish individual lines, rather one is pri- 
marily fitting the effects of a blend of lines (mostly iron L 
lines) at E ~ 1 keV. Where the emission actually consists of 
a mixture of gas phases spanning a range of temperatures, 
a single temperature fit can result in an underestimate of 
the metallicity. Similar low abundances have been seen in 
diffuse emission from a number of normal spiral galaxies 
(Read 1994), and in the case of elliptical galaxies, it has been 
reliably established by much higher spectral resolution stud- 
ies with the ASCA observatory, that the abundance of the 
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Figure 9. Surface brightness profile of the diffuse 0.1-1.9keV 
flux (Fig. 5) about the centre of the system. The enhancement at 
r ~ 25 kpc is due to the X-ray extensions. 



gas in hot X-ray emitting halos is ,$0.5 solar in many cases 
(Loewenstein et al. 1994; Matsushita et al. 1994). If the low 
abundance inferred for the gas in the Antennae is correct, 
then the X-ray emitting gas cannot represent the supernova 
ejecta, but must be derived from relatively uncontaminated 
gas from the galactic halos, or almost primordial material in 
the galaxies' environment. 

3.10 Overall properties and history 



The Antennae has a total 0.1-2.0 keV luminosity (leaving 
the system) of Lx = 1-1 x 10 41 erg s -1 , which is high com- 
pared to that of typical large spiral galaxies (Fabbiano 1989; 
David, Jones & Forman 1992; Read 1994). However a strong 
correlation is observed (Fabbiano, Gioia & Trinchieri 1988) 
between the X-ray and far-infrared luminosities among nor- 
mal and starburst galaxies, and the X-ray luminosity of the 
Antennaeis actually quite typical of high flux IRAS galaxies 
(Boiler et al. 1992). 

David et al. (1992) derive a relationship Lx °c 
-^Fifl" 1 01 (90% confidence limits) from Einstein data, 
though this is based upon rough approximations for the X- 
ray spectra (each spectrum is assumed to be a 5 keV ther- 
mal plasma absorbed by a column of 3 x 10 20 cm -2 ). This 
approximation may be appropriate for low activity spirals, 
which are dominated by emission from X-ray binaries, but 
in the case of the Antennae, which has a much softer spec- 
trum, it leads to a serious overestimate in its luminosity - 
David et al. (1992) quote an (0.5-4.5 keV) X-ray luminosity 
of 6.3 x 10 41 erg s" 1 , whilst our ROSAT results lead to a 
value of L x (0.5 - 4.5 keV) = 0.9 x 10 41 erg s" 1 . The lat- 
ter value compares well to the expected value from their 
X-ray/FIR correlation of 1.3 x 10 41 erg s" 1 . 



The conclusion that the Antennae actually has a normal 
X-ray luminosity for its activity, is confirmed by the fact that 
it falls close to the Lx-Lfir relationship derived for high 
FIR flux galaxies by Boiler et al. (1992), and to that derived 
with ROSAT iov a sample of nearby normal spiral galaxies 
by Read (1994). 

What of the way in which this X-ray emission is dis- 
tributed? The soft X-ray luminosity can be broken down 
into the following components: the discrete features A-G 
contribute 55%, the diffuse disc emission 41% (38% for the 
soft component attributed to hot gas, and 3% to the hard 
component from unresolved sources) and finally the two dif- 
fuse X-ray extensions account for 4% of the total. Sepa- 
ration of emission into discrete and diffuse components is 
sensitive to spatial resolution, and hence to source distance, 
and will also vary according to the analysis techniques em- 
ployed. It appears from a variety of studies with ROSAT 
(Walterbos et al. 1993; Bregman & Pildis 1994; Pietsch et 
al. 1994; Read 1994; Heckman 1993; Vogler, Pietsch & Ka- 
habka 1994; Schlegel 1994) that the fraction of genuinely 
diffuse emission (as distinct from harder emission from un- 
resolved sources) varies quite widely in spiral galaxies, cov- 
ering a range ~ 10 — 80%, but that a fraction ~ 50% is 
common, particularly in galaxies experiencing high star for- 
mation rates. The temperature of such diffuse emission is 
typically confined to the rather narrow band 0.3 — 0.6 keV. 
Hence on all these counts, the X-ray properties of the Anten- 
nae seem unremarkable. The only property which is unusual, 
is the morphology of the diffuse emission; combining a hot 
gas halo surrounding the galactic discs, with what appears 
to be a distorted but well-collimated outflow. 

On the basis of the considerations above, one can con- 
struct a tentative picture for the evolution of the X-ray 
properties of the system as follows. It seems likely that the 
X-ray extensions correspond to a wind driven from one or 
both galaxies at the time, ~ 1 — 2 x 10 8 years ago, when 
intense star formation was first triggered in their nuclei by 
the interaction. The bubble of very hot gas generated by the 
nuclear starburst will have broken out of the galactic disc 
perpendicular to the plane in the usual way, and once such 
a channel has been established, the wind should continue to 
flow so long as gas inflow within the disc continues to fuel 
further nuclear star formation. However the orientation of 
the outflows has been distorted from the normal oppositely 
directed bipolar form by the strong dynamical interaction 
of the two galaxies. This interaction has also been respon- 
sible for triggering distributed active star formation within 
the galactic discs, seen in the widespread Ha knots and also 
in the presence of luminous X-ray sources within the discs 
arising from emission from superbubbles and massive X-ray 
binaries. This star forming activity generates localised out- 
flows of hot gas into the galactic corona through 'chimneys' 
(Norman & Ikeuchi 1989), and this, together with further 
mixing and heating due to collisions between the gas associ- 
ated with each galaxy, has generated the hot corona which 
envelopes the system. 



4 CONCLUSIONS 

The observations reported here represent the first X-ray 
spectral study of a major galaxy merger in progress with 
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sufficient spatial resolution to resolve the various contrib- 
utors to the emission. The Antennae has a high X-ray lu- 
minosity compared to normal spirals, but falls well on the 
Lx-Lfir relationship established for late-type galaxies. A 
number of discrete components to the X-ray emission have 
been identified, and some of these have counterparts at other 
wavelengths. X-ray sources are seen near the two galactic nu- 
clei, but appear to derive from hot extranuclear gas rather 
than from the central starbursts themselves. In addition, X- 
ray emission is seen from two or three of the prominent Ha 
knots within the galactic discs. 

Nearly 50% of the flux associated with the system is 
diffuse, and its spectral properties suggests that the majority 
of this arises from hot gas rather than unresolved sources. 
The temperature of this gas is typically 4 x I0 6 K, and the 
total mass involved may be over I0 9 Mq, depending upon 
its filling factor. The bulk of the diffuse emission covers the 
optical discs of the two galaxies, and appears to envelope 
the whole optical system, apart from the tidal tails. This 
halo is cooler and less extensive than those typically seen 
around bright ellipticals but is substantially brighter than 
halos seen in normal spiral galaxies. An enhancement in the 
surface brightness and temperature of the diffuse emission 
is seen in the region where the two discs are believed to be 
colliding. This may result from direct heating of gas by the 
collision, or from star formation triggered in the vicinity. 

Two elongated structures extend to the NW and SW 
from the galactic discs to radii of ~ 30 kpc. It seems most 
likely that these indicate the presence of galactic winds, with 
a morphology which has been disturbed from the normal 
bipolar form by the rapid dynamical evolution of the system. 
They appear to contain approximately 10% of the total hot 
gas present, and could represent a significant loss of mass 
and energy from the system, though not of heavy elements 
if the low metallicities we infer are correct. 

Galaxy mergers such as we are witnessing in the An- 
tennae may have been much more common in the past (e.g. 
Larson (1990) and references therein) and could even be re- 
sponsible for a large fraction of elliptical galaxies. In this 
context it is interesting that we see what may be a hot gas 
halo in the making, and also substantial mass loss in a wind. 
Similar observations of systems in more advanced merger 
stages are required to chart the evolution of the hot gas 
through the transformation from spiral to elliptical galax- 
ies. 
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